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INTRODUCTION 


Both  finite  element  and  modified  Mapping  collocation  Mathoda  have  bean 
used  to  obtain  accurate  stress  Intensity  (K)  solutions  for  pressurized 
autofrettaged  thick  cylinders  Mlth  radial  cracks  (refs  1,2).  Use  of  the 
weight  function  Method  has  extended  two-dimensional  K  solutions  to  More 
refined  Material  Models  including  reverse  yielding  caused  by  the  Bauschinger 
effect  (ref  3).  Several  papers  have  used  stress  intensity  factors  to  estimate 
fatigue  lives  of  cannon  tubes  (refs  4-6).  The  calculations  underestimate  the 
measured  lives  for  pressurized  cylinders,  while  they  overestimate  the  experi¬ 
mental  results  for  pressurized  and  autofrettaged  cylinders  (ref  5).  The 
disagreement  between  measured  and  calculated  lives  diminishes  in  Reference  6 
by  introducing  a  fraction  of  the  negative  portion  of  K  values  as  a  part  of  the 
K  range. 

The  shallow  crack  approximations  for  K  solutions  were  used  and  a  linear 
approximation  for  Bauschinger  effect  on  residual  hoop  stress  was  assumed  in 
References  5  and  6.  The  accurate  two-dimensional  K  solutions  affected  by  a 
significant  Bauschinger  effect  using  elastic-plastic  analysis  were  used  in 
Reference  7  to  indicate  the  drastic  effect  of  reverse  yielding  on  the  life 
prediction.  Neither  shape  factors  nor  a  fraction  of  negative  K  were  con¬ 
sidered  in  Reference  7. 

In  this  report  the  life  prediction  formula  similar  to  that  used  in 
Reference  7  is  employed  to  check  the  calculated  lives  with  experimental 
results  of  Throop  (ref  8)  and  Throop  and  Fujczak  (ref  9).  The  stress  inten¬ 
sity  factors  for  surface  cracks  of  elliptical  shape  are  approximated  by  the 

References  are  listed  at  the  end  of  this  report. 
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two-dimensional  stress  intensity  factors  obtained  In  Reference  3  Multiplied  by 
respective  shape  factors  for  pressure  and  for  residual  stress  given  In 
Reference  5.  The  fraction  of  negative  K  Included  In  K  range  varies  froe  one 
at  the  notch  boundary  to  zero  at  a  crack  depth  far  away  from  the  notch.  The 
variation  of  this  fraction  Is  assumed  to  be  1/r*  where  r  Is  the  distance 
between  the  fatigue  crack  front  and  the  notch  front.  Another  Modification  is 
to  use  an  Initial  crack  depth  Much  deeper  than  the  notch  depth.  This  is  to 
avoid  large  cycles  required  in  experlMents  to  Initiate  a  single  continuous 
crack  front  along  the  notch  boundary  (it  is  considered  likely  (ref  5)  that 
Multiple,  small,  semi-elliptical  cracks  are  initiated  along  the  notch  boundary 
prior  to  their  link  together  to  form  a  single  crack).  With  these  modifica¬ 
tions  the  calculated  lives  agree  reasonably  well  with  experimental  results  for 
all  three  crack  shapes:  long  curves,  semi-el  1 iptical ,  and  semi-circular  used 
in  References  8  and  9. 

FATZQUE  LIFE  PREDICTION 

Integration  of  the  Paris  formula 

da 

--  -  C(AK)»  (1) 

dN 


for  fatigue  growth  rate  of  a  crack  subjected  to  cyclic  loading  is  usually  used 
to  determine  the  number  of  fatigue  cycles  required  tc  grow  a  crack  from  an 
initial  depth  a-j  to  a  final  depth  af. 


I 

N  «  Nf  -  Ni 


,af  d® 

*i  C(AK)m 


(2) 


where  C  and  m  are  material  constants  and  AK  is  the  range  of  stress  intensity 
defined  by 

*  *wax  "  *min  (3) 


*max  an<*  K*in  are  maximum  and  minimum  values  of  K  In  a  loading  cycle.  Assume 
that  a  crack  face  Is  a  geometric  plana  and  there  is  no  possibility  of  Inter¬ 
penetration  under  compression.  This  leads  to  a  conclusion  that  cannot 
be  a  negative  value.  In  the  case  of  repeated  firing  of  cannon  tubes 

«min  -  0  (4) 

is  used  for  both  autofrettaged  and  nonautofrettaged  tubes.  If  Kp  and  Kr 
denote  mode  I  K  values  corresponding  to  an  Internal  pressure  and  a  residual 
hoop  stress  respectively,  then 

*max  *  *p  (5a) 

for  nonautofrettaged  cylinders,  and 

*max  “  *p  *  *R  (5b) 

for  autofrettaged  cylinders.  Kp  and  Kr  are  usually  expressed  in  a  dimen¬ 
sionless  form  denoted  by  fKp  and  f*R,  respectively, 

*p  *R 

fKp  *  ZZ  *  fKR  ■  — ”  (6) 

prna  0o™a 

where  oQ  is  the  yield  stress,  p  is  the  Internal  pressure  applied  to  a  tube, 
and  p  is  related  to  aQ  by  a  load  factor  f|_  ■  o0/p.  By  virtue  of  Eqs.  (3), 

(4),  and  (5b) 

fKp  _ 

AK  «  (-—  +  (?) 

‘  L 

For  a  small  crack  growth  from  aj  to  aj+j,  f|(p  and  f*R  are  assumed  to  be 
constants  which  are  taken  as  the  mean  values 


*Kp  *  *(fKp<aj)  +  fKp(aj+l)) 
*KR  *  *(fKR(aj)  ♦  fKR(*j+l)) 


(8) 


The  integration  of  Eq.  (2)  becomes 


»  •  -  I  2-*(;!E  ♦  -  aj,,!'*/!)  (») 

Co  Co  j.1  fL 

where 

C°  *  c?i“5I 

with  t  denoting  the  Mali  thickness  and  a  ■  a/t.  The  fraction  a  is  used  since 
the  crack  depth  is  usually  expressed  as  a  fraction  of  Mall  thickness  t. 

Substituting  from  calculated  values  of  f^p  and  fgp  (ref  3)  corresponding  to 
various  material  models  in  Eqs.  (8)  and  (9),  Me  obtain  the  fatigue  crack 
groMth  (a  versus  N/Ca)  graph  shoMn  in  Figure  1  for  a  single,  two-dimensional 
straight-fronted  through  crack  in  a  cylinder  of  diameter  ratio  tMo.  In  this 
figure.,  f  is  the  Bauschinger  effect  factor.  The  dotted  lines  are  for 
idealized  material  Mithout  Bauschinger  effect  (f  -  1),  while  the  dashed  lines 
are  for  f  ■  0.38  (100  percent  overstrain)  and  f  ■  0.44  (60  percent 
overstrain),  respectively.  The  dashed  lines  with  m'  •  0  correspond  to 
elastic-perfectly  plastic  behavior  during  reverse  yielding.  The  dashed  lines 
with  m'  »  0.3  indicate  the  difference  in  predicted  cycles  when  strain¬ 
hardening  is  considered  in  reverse  yielding.  The  graph  shows  the  significant 
difference  of  autofrettage  effect  of  various  material  models  on  fatigue 
cycles.  Such  a  drastic  difference  is  not  supported  by  experimental  results. 
Corrective  parameters  must  be  introduced  to  correlate  the  calculated  fatigue 
cycles  with  observed  ones  in  laboratory  tests. 

SHAPE  FACTORS 

The  ratio  of  stress  intensity  factor  for  a  surface  crack  to  that  of  a 
two-dimensional  through  crack  is  called  the  shape  factor.  The  stress  inten¬ 
sity  factor  varies  along  the  crack  front  of  a  surface  crack.  It  changes  with 
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crack  shape  and  under  different  loadings.  If  the  variation  of  stress  Inten¬ 
sity  along  the  crack  front  Is  Important ,  the  three-dimensional  X  solutions  for 
the  surface  crack  should  be  obtained.  If  an  estimate  of  K  at  a  point,  say  the 
deepest  point  of  the  surface  crack.  Is  needed,  then  a  reasonable  estimate  of 
shape  factor  Is  useful.  An  extensive  study  of  shape  factors  has  been 
published  by  Hetman  and  Raju  (ref  10)  for  semi-elliptical  cracks  In  a  flat 
plate  under  tension  or  bending.  There  Is  no  comparable  study  for  such  a  crack 
In  a  thick  cylinder.  Parker  et  al  (ref  S)  obtained  estimates  of  shape  factors 
for  semi-elliptical  cracks  of  various  aspect  ratios  In  a  pressurized  and  auto- 
frettaged  cylinder  from  judicious  use  of  results  reported  In  Reference  10. 

More  accurate  three-dimensional  K  solutions  should  be  obtained  to  check  these 
estimates.  Before  more  accurate  and  reliable  shape  factors  become  available, 
values  close  to  those  given  by  (a)  and  (c)  of  Figure  7  In  Reference  5  are  used 
for  fSp  and  f$R.  respectively,  in  this  study.  Multiplying  fKp  and  fKR  with 
their  shape  factors,  Eq.  (7)  becomes 

AK  ■  (-3J-SB  +  fSRfKR)«70/ira  (11) 

fl 

The  following  mean  values  are  used  in  Eq.  (11)  for  a  small  crack  growth  from 
•j  t0  aj+l* 

^Sp^Kp  *  J*tfsp(*j)^Kp(®j)  +  ^Sp(aj+l^Kp(*j+:i)3 

_  (12) 

fSRfKR  ■  *£fSR<aj)fKR(aj)  ♦  fSR<aj+l )fKR(aj*l > 3 

NEGATIVE  K  FACTOR  FN 

For  a  fatigue  crack  which  is  considered  as  a  geometrical  plane  with  no 
thickness,  Kn^n  ■  0  is  used  in  Eq.  (3)  to  obtain  K  range.  For  a  notch  with 
finite  thickness,  the  upper  and  lower  planes  of  the  notch  may  have  normal 
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displacement  In  both  directions.  The  argument  used  to  limit  a  0  Is  not 
valid  at  the  notch  front.  In  fact,  the  full  negative  K  should  be  used  for  a 
crack  starting  at  the  notch.  At  some  depth,  the  notch  effect  nay  become 
negligibly  small,  and  KM^n  ■  0  nay  again  be  used  In  Eq.  (3)  for  AK.  Kendall 
has  argued  that  a  portion  of  the  negative  K  must  be  included  In  calculating  K 
range  (ref  6).  He  introduced  a  constant  fraction  which  is  multiplied  by  the 
negative  K  (Kp)  coresponding  to  compressive  rnsidual  stress  due  to  autofret- 
tage  to  give  the  Km^n  in  AK.  From  our  hypothesis,  the  fraction  varies  with 
the  depth  of  the  fatigue  crack.  As  a  first  approximation,  the  fraction  fp  is 
taken  as  (1  ♦  r/rn)“«  where  rn  is  the  depth  of  the  notch  end  r  is  the  depth  of 
the  crack  measured  from  the  notch  front.  At  the  notch  front  r  ■  0  and  fN  a  l, 
the  full  negative  K  is  taken  as  Km-jn.  When  r  ■  rn  the  fatigue  crack  grows  to 
a  depth  equal  to  the  notch  depth;  this  approximation  gives  fp  *  fc .  It  is  also 
assumed  that  fp  ■  0  when  r  >  2rn.  A  linear  variation  of  fp  was  another 
approximation  examined  which  underestimated  the  measured  fatigue  lives. 
Incorporating  fp,  Eq.  (11)  becomes 


'Kp  , — 

■  [fsp  I*’  ♦  U-fN)fSRfKRW** 

'  L 


Equation  (9)  can  be  expressed  explicitly 

N  "  *  fKpdj.l)  -  -  1-./2  l-m/2.  .... 

Z  *  l  t . . . . +  fKR(«j)  +  fKR(®j+l)l  -  ci  +  1  )  H4) 

C0  j»i  fL  J  J  i 

where  abbreviations  f«p,  f are 

fKp(«j)  ■  fsp(aj)fKp(aj) 

(15) 

fKR(aj)  «  [1  -  fN(0j>]fSR(°j)fKR(<*j) 


LABORATORY  SPECIMENS  AND  MEASURED  LIVES 


Throop  (ref  8)  and  Throop  and  Fujczak  (ref  9)  obtained  a  series  of 
experimental  results  which  relate  fatigue  crack  lives  with  crack  shape  and 
extent  of  autofrettage  for  thick-wall  cylinders.  The  cylinders  were  0.76  m 
(30  inches)  in  length,  180  mm  (7.1  inches)  bore  diameter,  360  mm  (14.25 
inches)  outside  diameter  and  were  fatigue  cracked  from  longitudinal  internal 
notches.  Three  initial  notch  geometries  were  used;  semi-circular,  100  mm 
(4-inch)  and  500  mm  (20-inch)  long  notches  produced  by  electrical  discharging 
machining.  They  were  6.4  mm  (%-inch)  deep  by  0.76  mm  (0.03-inch)  wide,  the 
semi-circular  notch  being  13  mm  (3<-inch)  diameter  half-penny  shape.  Fatigue 
cracks  grown  from  the  initial  notches  were  monitored  periodically  for  depth 
and  shape  with  ultrasonics  as  the  cylinder  was  repeatedly  pressurized  to  330 
MPa  (48  Ksi).  The  cylinder  material  was  ASTM  A723  forged  steel,  with  yield 
strength  of  1175  MPa,  -40#C  Charpy  impact  energy  of  34  J,  reduction  in  area  of 
50  percent.  A  schematic  diagram  of  a  typical  cylinder  with  a  simple  initial 
notch  and  the  growth  of  a  500  mm  (20-inch)  long  notch  from  6.4  mm  (Vinch) 
initial  depth  by  repeated  pressurization  is  shown  in  Figure  2.  The  measured 
crack  depth  versus  corresponding  number  of  fatigue  loadings  is  shown  in  Figure 
3  for  a  single-notched  cylinder  with  0,  30,  and  60  percent  overstrains  for 
three-notch  geometries. 

Since  the  fatigue  crack  is  not  likely  to  start  from  the  notch  imme¬ 
diately,  an  adjustment  of  experimental  results  is  made  by  subtracting  the 
number  of  cycles  required  to  grow  from  initial  notch  depth  to  an  initial  crack 
depth  (a.j)  from  the  number  of  cycles  to  grow  from  initial  notch  depth  to  a 
final  crack  atpth  (af).  The  initial  crack  depth,  which  should  be  reasonably 
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larger  than  the  initial  notch  depth  (6.4  mm  in  this  experimental  wtudy),  is 
arbitrarily  taken  as  a-j  ■  O.lt  ■  9  mm.  The  original  experimental  data  are 
available  only  in  graphs.  To  improve  the  accuracy  of  readings  from  graphs, 
the  graphs  were  first  enlarged  and  the  average  was  used  from  values  obtained 
from  different  graphs  published  in  different  papers  (refs  4,5,9)  for  the  same 
set  of  experimental  data.  The  adjusted  experimental  results  are  shown  by  dots 
in  Figures  4,  5,  6,  and  7. 

PREDICTED  LIVES  FOR  THE  EXPERIMENTS 

The  material  constants  for  the  steel  used  in  the  experiments  are  m  =  3 
and  C  *  6.52xl0-1 *  for  crack  growth  in  meters  per  cycle  and  AK  in  MPa/meter 
(or  C  *  3.4xl0-10  for  crack  growth  in  inches  per  cycle  and  AK  in  Ksi/Tnch). 
Using  f^  ■  3.55  and  values  of  fgp  and  fgp  in  Table  I,  and  assuming  = 

(1  +  r/rn)~*  for  0  <  r  <  2rn  and  fp  «  0  for  r  >  2rn,  Eqs.  (14)  and  (15)  can  be 
evaluated  with  known  discrete  values  of  fKp  and  f kr.  Values  of  fgp  and  fgR, 
given  in  Table  I,  are  obtained  from  Reference  5  for  different  crack 
geometries.  Discrete  values  of  f«p  and  f«R  can  be  computed  by  the  method 
described  in  Reference  3.  The  calculated  lives  and  measured  lives  are  plotted 
in  Figures  4,  5,  and  6  for  three  crack  geometries,  respectively.  For 
nonautofrettaged  cylinders  (zero  percent  overstrain),  there  is  only  one  set  of 
calculated  results,  shown  in  a  solid  line,  for  each  crack  configuration.  For 
60  percent  overstrain,  solid  lines  are  for  idealized  material  without  con¬ 
sidering  the  Bauschinger  effect,  while  two  dashed  lines  in  each  figure  are 
predicted  lives  with  some  reverse  yielding.  The  two  dashed  lines  differ  in  m' 
values,  m'  ■  0  or  0.3,  where  m'E  is  the  slope  of  strain-hardening  during 
reverse  yielding. 
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When  m'  ■  0,  the  material  behaves  like  elastic-perfectly  plastic  in 
reverse  yielding.  The  magnitude  of  the  compressive,  residual  stress  in  the 
tangential  direction  near  the  bore  is  smaller  than  that  in  a  strain-hardening 
material  (ref  3).  The  strain-hardening  reduces  the  adverse  effect  of  reverse 
yielding  due  to  the  Bauschinger  effect.  Figure  1.  Figures  4  through  6  show  an 
overall  effect  of  all  factors.  Figure  7  shows  the  effect  of  each  factor  on 
the  semi-elliptical  surface  crack  in  a  pressurized  cylinder  with  60  percent 
overstrain.  The  dotted  curve,  curve  1,  gives  the  calculated  lives  for  an 
idealized  material  with  no  Bauschinger  effect.  No  correction  factors  are 
used.  The  stress  intensity  range  is  based  on  two-dimensional  stress  intensity 
factors  and  Km^n  =  o.  If  the  change  in  residual  hoop  stress  due  to  the 
Bauschinger  effect  (Bauschinger  effect  factor  f  =  0.44)  is  considered,  the 
calculated  lives  are  shown  as  a  dashed  curve,  curve  2,  where  m'  =  0.3  is  used. 
The  large  difference  between  curves  1  and  2  shows  the  significant  effect  of 
reverse  yielding  on  fatigue  lives.  If  shape  factors  fgp  and  fgR  are  con¬ 
sidered  in  addition  to  the  Bauschinger  effect,  the  predicted  lives  will  be 
changed  from  curve  2  to  curve  3.  The  shape  factors  are  to  increase  the 
fatigue  lives.  Finally,  if  the  correction  factor  f^  is  also  taken  into  con- 
sidaration,  the  predicted  lives  are  shown  as  the  solid  line,  curve  4,  which  is 
close  to  the  experimental  results  shown  by  dots  in  Figure  7.  Similar  graphs, 
obtained  for  the  long  curved  crack  and  the  semi-circular  crack,  are  omitted 
since  they  indicate  similar  effects  of  each  correction  factor. 
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TABLE  I.  VALUES  OF  fsp  AND  fSR  FOR  DIFFERENT  CRACK  GEOMETRIES 


o  *  a/t 

20- Inch 

Lona  Notch 

4 

Lon 

-Inch 
a  Notch 

Semi-Circular 

Notch 

fsP 

fSR 

fsp 

fSR 

fSp 

fSR 

0.1 

0.95 

0.95 

0.72 

0.69 

0.57 

0.58 

0.2 

0.90 

0.90 

0.64 

0.60 

0.53 

0.52 

0.3 

0.85 

0.85 

0.61 

0.55 

0.51 

0.42 

0.4 

0.85 

0.80 

0.57 

0.47 

0.48 

0.34 

0.5 

0.85 

0.75 

0.56 

0.38 

0.44 

0.24 

CONCLUSIONS 

The  fatigue  life  of  a  thick-walled  cylinder  can  be  predicted  reasonably 
well  by  integration  of  the  Paris  formula  for  growth  rate  of  a  fatigue  crack 
under  cyclic  loading,  if  a  modified  stress  intensity  range  is  used.  The 
stress  intensity  range  is  obtained  by  multiplying  two-dimensional  stress 
intensity  factors  by  proper  shape  factors  and  some  negative  K  factors.  More 
systematic  and  controlled  experiments  are  required  to  check  the  idea  proposed 
in  this  report.  Three-dimensional  K  solutions  for  semi-elliptical  surface 
cracks  in  thick  cylinders  with  various  residual  stresses  are  needed  to  esti¬ 
mate  the  proper  shape  factors.  Special  experiments  should  be  performed  to 
verify  the  concept  of  negative  K  factors  and  to  determine  the  variation  of  f^ 
in  terms  of  crack  depth. 
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Figure  1.  Effect  of  reverse  yielding  on  fatigue  crack  growth  shown  in  a 
(a  versus  N/C0)  plot  for  a  2.0  diameter  ratio  cylinder  with  a 
single,  straight-fronted  through  crack. 


N,  (IOOOCYCLES)  N.  (IOOO  CYCLES ) 

(a)  LONG  CURVED  CRACK  (b)  SEMI-ELLIPTICAL  CRACK 


(c)  SEMI-CIRCULAR  CRACK 

Figure  3.  Creek  depth  plotted  versus  number  of  fatigue  loadings 
experimental  results  (ref  9). 
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